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Nucleation, commonly associated with discontinuous transformations between metastable and stable 
phases, is crucial in fields as diverse as atmospheric science and nanoscale electronics. Traditionally, it is 
considered a microscopic process (at most nano-meter), implying the formation of a microscopic nucleus of 
the stable phase. Here we show for the first time, that considering long-range interactions mediated by 
elastic distortions, nucleation can be a macroscopic process, with the size of the critical nucleus proportional 
to the total system size. This provides a new concept of "macroscopic barrier- crossing nucleation". We 
demonstrate the effect in molecular dynamics simulations of a model spin-crossover system with two 
molecular states of different sizes, causing elastic distortions. 

Nucleation is a barrier- crossing process 1 " 3 , in which a metastable phase decays via a critical nucleus for 
which the increase in surface free energy is compensated by the bulk energy decrease 1 " 10 . If the cluster 
becomes bigger than the critical size, it grows, while if smaller, it shrinks. The size of a critical nucleus is 
determined by microscopic competition between the surface and bulk free energies of a microscopic cluster, and 
thus the size of the critical nucleus is microscopic 3 ' 7 " 9 (see Fig. la and b). To be precise, this situation is realized in 
short-range interaction systems, where separation of the energy between the bulk and surface is allowed. 
Nucleation theories 1 " 410 have been based on this idea, and so far only microscopic nucleation is known. 

Consider a typical short-range interaction system at low temperatures: the d- dimensional Ising model 11 
defined by the Hamiltonian, H=JJ2 i j(Ji(Jj — hJ2i (7 i> where a = ±1 (up and down spins). The free-energy 
barrier for a droplet with radius r is AE = — C h hr d + C s Jr d ~\ Here C h r d is proportional to the volume of the droplet 
and C s r d ~ l is proportional to the area of the phase boundary. As depicted in Fig. lb, the radius of the critical 
droplet is given when the droplet has the maximum excess free energy as r c = c^jh ■ The critical radius r c is 
independent of the system size. Figure lc shows an example of nucleation in a circular system (open boundary 
conditions, OBC) for the Ising model (d = 2). Nucleation takes place both in the bulk (inside) and at the 
boundary. Nucleation at the boundary is energetically more favorable, but when the system becomes larger 
(the bulk-to-boundary ratio becomes large), nucleation in the bulk becomes dominant (Supplementary 1). 

However, when the interaction is of long range, the nature of nucleation is different. Because it has been pointed 
out that an elastic interaction due to lattice distortion causes an effective long-range interaction 12 , the nucleation 
process in systems with elastic interactions (e.g., spin-crossover systems 13 " 16 , martensitic systems 17 " 20 and Jahn- 
Teller systems 21 " 24 ) should be investigated. In this work we present properties of the nucleation in long-range 
elastic interaction systems with OBC. 

In molecular crystals, e.g., transition-metal complexes, a molecule often displays bistability in both its electro- 
nic state and molecular size (structure). External stimuli, e.g., change of temperature, pressure, photoirradiation, 
etc. change the molecular size. The distortion caused by the change of size induces an elastic interaction, which 
acts as an effective long-range interaction 12 ' 25,26 . Spin -crossover (SC) compounds are a typical example of the 
above situation (see Fig.2 a-c), where the low-spin (LS) and high-spin (HS) states are separated by an energy 
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barrier, and the LS molecule is smaller than the HS one. Indeed, SC 
systems show a wide variety of phase transitions under external 
stimuli 13 " 16 . 

In the present work we study the nucleation dynamics of circular 
(d = 2) crystals of a long-range elastic interaction system and show 
that the nucleation is a barrier-crossing process. However, the size of 
the critical nucleus (r c ) is proportional to the system size (R). Thus, a 
macroscopic nucleation mechanism is realized, which is qualitatively 
different from previously known nucleation mechanisms. 



C 




Figure 1 | Nucleation and domain formation for short-range interaction 
systems, (a) A schematic example of a droplet of a short-range magnetic 
interaction system. The blue circle with radius r is a cluster of down- spin 
molecules in an up-spin phase (red part), (b) Microscopic competition of 
the surface and bulk free energies of a droplet. The surface free energy is an 
increasing function of the radius of the droplet ( r) (upper broken line) and 
the bulk free energy is a decreasing function of r (lower broken line). The 
solid line is the sum of these energies (AE). The critical nucleus with the 
radius r c gives the maximum of the total free energy, (c) Snapshots of 
nucleation for a short-range interaction model (Ising model). Red and blue 
denote up and down spins, respectively. 
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Results 

We adopt the following Hamiltonian for the model 27 , 

N ^2 N N p2 

n 0 =y^ + y y; ntra (^-) + Y — 
+ E y T r ( x ^>^)- 

The first and second terms describe the motion of the intramolecular 
mode of the zth molecule. The radius of the molecule is r b and the 
conjugate momentum is p t . The mass for the motion is m 
(Supplementary 2). The intramolecular potential energy Vf Am {x i ) 
is shown by the solid curve in Fig. 2a, where x { = r,- — r LS . Here r LS 
is the ideal radius of the LS molecule, and that of the HS molecule is 
Ths = r LS + Ar. Vj ntm (xi) provides a symmetric vibration mode and 
induces changes of the molecular size. The third and fourth terms 
describe the ceneter-of-mass motion of the molecules (X b P { ) with 
mass M. The intermolecular interaction V inter (see Methods) for 
nearest and next-nearest neighbors is a function not only of the 
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Figure 2 | Elastic interaction system and nucleation features, (a) Intramolecular potential energy V(x) shown by the solid (blue) curve. The dotted 
curves are low spin (LS) and high spin (HS) potential energies without quantum mixing. The curvature for the LS state is 4 times larger than for the HS 
state 13 in this work but other choices for the ratio of the curvatures, for example equal curvatures, do not change the essence of the results. The energy unit 
is 100-300 K for SC compounds 13 (Supplementary 2), and it is also the unit of the temperature T. (b) LS molecule (blue) and HS molecule (red). The HS 
molecule is larger in size, (c) Distortion due to the difference of the molecular sizes, (d, e) Snapshots of the configuration during relaxation from the HS 
phase for (d) 2R = 200 and (e) 2R = 100. The value of/ HS is 0.95 for dl and el, 0.90 for d2 and e2, 0.71 for d3 and e3, 0.50 for d4 and e4, 0.26 for d5 and e5, 
and 0.15 for d6 and e6. 
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coordinates X { and Xp but also of the molecular radii r { and r ; . 
Although the interaction terms in this Hamiltonian appear only 
to be of short range, elastic interactions mediate the effect of the local 
lattice distortions over long distances 12 . We study relaxation and 
nucleation in this model at low temperatures by using a molecular 
dynamics method 28 . The critical temperature of the first- order trans- 
ition of the model is T c ~ 0.9, below (above) which the LS (HS) state is 
the equilibrium state 25 . 

We observe relaxation from the metastable HS state to the LS state 
at a low temperature (T = 0.2) in approximately circular crystals 
on a square lattice. This temperature is much lower than the critical 
temperature T c . For the initial states of the relaxation (the metastable 
HS phase), we gave a set of velocities to all molecules according to 
the Maxwell- Boltzmann distribution by using a random number 
sequence. 

Snapshots of the configuration during the course of a relaxation 
event are depicted in Fig. 2d and Fig. 2e, where the diameters of the 
circular crystals are 200 and 100 particles (denoted as 2R = 200 and 
2R = 100), respectively. Figures 2dl and 2el show configurations 
when the HS fraction (f H s) 28 reaches the value / H s — 0.95 for 2R = 
200 and 2R = 100, respectively. Nucleation starts from one point 
along the circumference. The subsequent configurations are given in 
Figs 2d(2-6) and 2e(2-6) for 2R = 200 and 2R = 100, respectively. 
The corresponding values off H s are the same in both systems. 



As we show below, the configurations of Figs 2dl and 2el are those 
of the critical nucleus, and Figs 2d(2-6) and Figs 2e(2-6) correspond 
to deterministic growth of the LS droplet after the formation of the 
critical nucleus. It should be noted that the LS domain shapes can be 
well characterized by using the contact angle (wetting angle) 29 of nil. 
We checked the configurations for other relaxations from different 
initial conditions (different random number sequences for the 
molecular velocities) and found the same features of nucleation 
and growth. Here the size of the critical nucleus is found to be 
proportional to the system size. The shapes of the critical nuclei 
and also the following clusters are geometrically similar in systems 
of different sizes. This feature is qualitatively different from that of 
short-range interaction systems 30 , in which the critical droplet has a 
specific size independently of the system size. 

In Figs 3a and 3b the time dependence of/ H s is shown for sys- 
tems with 2R = 100 and 2R = 200, respectively. Because a single 
nucleation event dominates the process, the escape time from the 
metastable state is random and governed by a Poisson process. 
However, once nucleation starts, the process is almost deterministic. 
The crossover between the stochastic and deterministic regimes 
determines the critical nucleus size and a threshold value of / HS . 
These are typical characteristics of barrier- crossing dynamics. 

To capture this feature, we study the relaxation of/ H s a ft er passing 
a given value of /hs(=/hs)- The passing time t tr is defined as 
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Figure 3 | Relaxation processes from the metastable HS phase, (a, b) HS fraction (f HS ) versus time at T = 0.2 for (a) 2R = 100 and (b) 2R = 200. The time 
of the collapse of the metastable state is widely distributed, (c, d) The three panels of (c) show the time dependence of fn S > starting from the time t tr , for 
2R = 100 when fg s =0.98, 0.95, and 0.9 (from left to right). The three panels of (d) correspond to the case of 2R = 200. 
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fm = fns(ttr)- As mentioned above, the time evolutions after passing 
the threshold value, i.e., fns(t — t tr ) are expected to overlap. We 
plotted the data of/ H s(f — t tr ) for various trial values of/^ r s , and 
found /ns(^ tr ) —0.95 gives the threshold as depicted in Fig. 3c. The 
same value is observed in both systems with 2R = 100 (Fig. 3c) and 
2R = 200 (Fig. 3d), and we conclude that it is independent of the 
system size. 

Discussion 

To examine the features of the critical nucleus and check the size 
dependence, we analyze the total potential energy of the system 
(£ tot = £ V™ tra + Vf ex ) as a function of the relative size of the LS 
domain. It is considered that the entropy effect is small enough 
compared to the energy barrier during the relaxation at this low T. 
As a parameter to characterize the domain size, we define 0 as 
the central angle. With the contact angle of n/2, the domain region 
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Figure 4 | Barrier-crossing of macroscopic nucleation. (a) The excess 
energy density Ap as a function of 6 for the system sizes 2R = 100, 200, and 
300. The inset is the definition of 6 and domain region is given for 6 using 
the contact angle of n/2. For all R, the values of / HS = 0.98, 0.97, 0.95, 0.92, 
0.90, 0.86, and 0.83 are given at 6/n = 1/10, 1.7/10, 2.3/10, 2.9/10, 3.5/10, 
4.2/10, and 4.8/10 respectively. After 6/n = 4.8/10, Ap decreases 
monotonically until 6/n = 2 (LS phase), (b) The dependence of Ap on 1/R 
at 6 = 2.3ti/10 for 2R = 100, 200, 300, 400, and 500. Ap approaches a finite 
value as R approaches infinity. 



(lens-shaped part) is defined for any 9 (Fig. 4a), where the inter- 
face between the two phases is given by the circle of the radius 
r d = R tan (0/2), whose center is the crossing point of the two 
tangential lines. The value of £ tot for a given 0 is obtained as follows. 
In the circle of the HS phase, we replace HS molecules in the lens- 
shaped part subtended by 6 by LS molecules. Then we move all 
molecules slowly so as to reach the minimum total potential - 
energy state, and obtain the energy value of this stationary state. 
We define the energy density as p = %S where N is the number 
of molecules in the system and N ~ tlR 2 , and also the relative energy 
density: Ap = p — p e=0 as the difference between p of the stationary 
state and that of the complete HS phase (po=o)> 

We show Ap as a function of 0 for several system sizes (2R) in 
Fig. 4a. For small values of 0(< 7i/10), Ap is almost constant and then 
Ap increases with 0. In this region the cluster is expected to shrink in 
the relaxation process. Around 0 = 2.37T/10, Ap shows the maximum 
value and it decreases for larger 0. It should be noted that at this 
0( = 2.37i/ 1 0) f HS is equal to 0.95, which agrees with the threshold 
value off H s hi the analysis of the relaxation curves (Fig. 3c and d). 

When the droplet size exceeds the critical size, the domain 
expands. For different system sizes, this critical size of the droplet 
(ra) changes in proportion to the system syze (R). Namely, the critical 
angle exists, but not a specific critical size. This fact was demonstrated 
in Fig. 2d and Fig. 2e, i.e., the domain shape is almost the same for 
systems of different size. Thus we call this process "macroscopic 
nucleation", and we believe that it should hold even in the bulk 
(continuum) limit. We depict Ap at 0 = 23n/l0 (peak position) 
as a function of 1AR in Fig. 4b, and find the dependence: 
Ap = Ap 0 — £2sst e The value Ap 0 ~ 0.035 is considered the value of 
the bulk limit. 

The behavior of "macroscopic nucleation" is qualitatively differ- 
ent from that observed in short-range interaction systems. The bulk 
and surface contributions to the potential barrier cannot be distin- 
guished in this long-range interaction system, which is similar 
to interface energies of binary alloys due to elasticity 31 , and the 
elastic interactions suppress both bulk nucleation and multi-droplet 
nucleation at the boundary (Supplementary 1 and 3). 

In summary, we propose a new concept of "macroscopic nuc- 
leation" for systems with long-range interactions. The domain 
formation exhibits geometric similarity for circular crystals of 
any size. This means that the size of the critical nucleus is propor- 
tional to the system size and macroscopic nucleation is realized. 
Recognition of this mechanism should give important insights for 
all systems in which local structural changes cause a distortion of the 
lattice. In addition to the spin-crossover type systems considered 
here, the mechanism should hold for martensitic and Jahn-Teller 
systems, etc. 

Methods 

The role of the intermolecular potential is to release the local distortions due to the 
difference of the sizes of neighboring molecules. For this purpose, we 
adopt the following potential 27 , Vf tet (X ii X j> r h rj) =f(d ij -Ar), where 
f(u) =D(e fl ' (u ~ Mo) +e-k>-"°)). The variable u 0 is a constant such that/(w) has its 
minimum at u = 0 and dy = |JQ — Xj\ — (r; + rj). For nearest neighbors, Ar = 0, 
a' = 0.5 and b' = 1.0 are set, and the energy minimum is realized when the 
neighbors have the same size. For next-nearest neighbors, Ar = 2(y / 2 — l)r with 
r= (r LS + r H s)/2, a' = 0.1 and b' = 0.2 are set. This provides a small force 
sufficient to ensure the stability of the crystal structure (this is specific to coor- 
dination z = 4). The parameter D associated with the strength of the intermol- 
ecular interaction was set to D = 20, which is strong enough to cause a first- order 
phase transition. Here the type of the potential function is not so important, and 
the basic mechanism of macroscopic nucleation is universal for other types of 
intermolecular potentials (harmonic or anharmonic potentials). The other para- 
meters were set as r H s = 9, Ar = 1, and m = M = 1 (Supplementary 2). Molecular 
dynamics simulations were performed using a Nose-Hoover thermostat 28 . With this 
method, the timescale of the simulation is influenced by the thermostat parameters. 
Here, we used this effect to our advantage to perform the simulations for large 
systems in a computationally feasible time. Although we sacrifice the ability to 
measure nucleation times, which we expect to increase dramatically with system size, 
our method allows us to observe the scale-invariant spatial structure of the nucleation 
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process in systems of very different sizes, as shown in Figs 2d and 2e. The qualitative 
feature of nucleation in a stochastic Poisson process, followed by deterministic 
growth shown in Fig. 3, is also preserved. 
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